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Abstract 



We have studied charmless hadronic decays of B mesons into two-body fi- 
nal states with kaons and pions. We present preliminary results based on 
9.66 million BB pairs collected with the CLEO detector. We have made 
the first observation of the decay B vr+vr", with the branching frac- 
tion of Br{B vr+vr-) = (4.7111 ± 0.6) x 10"^. We have also observed 
for the first time the decay B K^tt^ with the branching fraction of 
Br{B K^TT^) = (14.8l|?l||) x lO'^, thus completing the set of four Ktt 
branching fraction measurements. We present improved measurements for the 
decays B — > K^ir^, B^ —>■ K^tt^, and B^ —)■ K^tt^. We use these and other 
charmless hadronic B decays to make a first determination of the value of the 
weak phase Avg{V*^) = 7 = 113°l23°- 
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I. INTRODUCTION 



The phenomenon of CP violation, so far observed only in the neutral kaon system, can 
be accommodated by a complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark- 
mixing matrix [|l|. Whether this phase is the correct, or only, source of CP violation awaits 
experimental confirmation. B meson decays, in particular charmless B meson decays, will 
play an important role in verifying this picture. 

The decay B —>■ vr+vr", dominated by the b —>■ u tree diagram (Fig. |I](a)), can be used to 
measure CP violation due to B^ — B^ mixing at both asymmetric B factories and hadron 
colliders. However, theoretical uncertainties due to the presence of the b ^ dg penguin 
diagram (Fig. |l](b)) make it difficult to extract the angle a of the unitarity triangle from 
B TT+TT^ alone. Additional measurements of B^ vr^vr^, B, B ^ tt^tt'^, and the use of 
isospin symmetry may resolve these uncertainties [Q]. 

3460997-007 




(c) (d) 

FIG. 1. The dominant decay processes are expected to be (a) external W-emission, (b) gluonic 
penguin, (c) internal W-emission, (d) external electroweak penguin. 

B —>■ Ktt decays are dominated by the b ^ sg gluonic penguin diagram, with additional 
contributions from b ^ u tree and color- allowed electroweak penguin (Fig. |l](d)) processes. 
Interference between the penguin and tree amplitudes can lead to direct CP violation, which 
would manifest itself as a rate asymmetry for decays of B and B mesons. Several methods 
of measuring or constraining 7, the phase of Vub- using only decay rates oi B ^ Kn, nn 
processes were also proposed @] 0- This is particularly important, as 7 is the least known 
parameter of the unitarity triangle and is likely to remain the most difficult to determine 
experimentally. This paper presents the first observation of the decays B vr+vr" and 
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B /^°7r°, improved measurements for B K'^tt^, B^ K^tt^, and B^ K^ti^ 
decays, and updated upper limits for B decays to vr^vr^, K^K~ , and K^K^. Average over 
charge conjugate decays is implied throughout this paper. 

II. DATA SET, DETECTOR, EVENT SELECTION 

The data used in this analysis was collected with the CLEO II and CLEO II. V detectors 
at the Cornell Electron Storage Ring (CESR). It consists of 9.1 fb~^ taken at the T(4S) 
(on-resonance) and 4.5 fb~^ taken below BB threshold. The below-threshold sample is used 
for continuum background studies. The on-resonance sample contains 9.66 million BB pairs. 
This is a 67% increase in the number of BB pairs over the previously presented analysis 0. 

CLEO II and CLEO II. V are general purpose solenoidal magnet detectors, described in 
detail elsewhere [^]. In CLEO II, the momenta of charged particles are measured in a tracking 
system consisting of a 6-layer straw tube chamber, a 10-layer precision drift chamber, and 
a 51-layer main drift chamber, all operating inside a 1.5 T superconducting solenoid. The 
main drift chamber also provides a measurement of the specific ionization loss, dE/dx, used 
for particle identification. For CLEO II. V the 6-layer straw tube chamber was replaced by 
a 3-layer double-sided silicon vertex detector, and the gas in the main drift chamber was 
changed from an argon-ethane to a helium-propane mixture. Photons are detected using 
7800-crystal CsI(Tl) electromagnetic calorimeter. Muons are identified using proportional 
counters placed at various depths in the steel return yoke of the magnet. 

Charged tracks are required to pass track quality cuts based on the average hit residual 
and the impact parameters in both the r — and r — z planes. Candidate Kg are selected from 
pairs of tracks forming well-measured displaced vertices. Furthermore, we require the Kg 
momentum vector to point back to the beam spot and the 7r"'"7r~ invariant mass to be within 
10 MeV, ~ 2.5 standard deviations (a), of the Kg mass. Isolated showers with energies 
greater than 40 MeV in the central region of the Csl calorimeter and greater than 50 MeV 
elsewhere, are defined to be photons. To reduce combinatoric backgrounds we require the 
lateral shapes of the showers to be consistent with those from photons. To suppress further 
low energy showers from charged particle interactions in the calorimeter we apply a shower 
energy dependent isolation cut. Pairs of photons with an invariant mass within 25 MeV 
(~ 2.5cr) of the nominal tt" mass are kinematically fitted with the mass constrained to the 
tt" mass. 

Charged particles [h^) are identified as kaons or pions using dE/dx. Electrons are re- 
jected based on dE/dx and the ratio of the track momentum to the associated shower energy 
in the Csl calorimeter. We reject muons by requiring that the tracks do not penetrate the 
steel absorber to a depth greater than seven nuclear interaction lengths. We have studied 
the dE/dx separation between kaons and pions for momenta p ~ 2.6 GeV/c in data using 
i5*^-tagged K~tt~^ decays; we find a separation of (1.7 ± 0.1) a for CLEO II and 

(2.0 ±0.1) a for CLEO ILV 
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III. ANALYSIS 



We calculate a beam-constrained B mass M = y E'^ — p^, where pb is the B candidate 
momentum and i?b is the beam energy. The resolution in M is dominated by the beam 
energy spread and ranges from 2.5 to 3.0 MeV/c^, where the larger resolution corresponds 
to decay modes with a 7r°. We define AE = Ei + E2 — Ey,, where Ei and E2 are the energies 
of the daughters of the B meson candidate. The resolution on AE is mode dependent. 
For final states without 7r°'s, the AE resolution for CLEO II(II.V) is ±25(20) MeV. In 
the B^ h^ir^ analysis the AE resolution is worse by approximately a factor of two and 
becomes asymmetric because of energy loss out of the back of the Csl crystals. The energy 
constraint also helps to distinguish between modes of the same topology. For example, AE 
for B K^TT~, calculated assuming B tt^tt^, has a distribution that is centered at 
—42 MeV, giving a separation of 1.7(2.1)cr between B — >■ K~^7r~ and B — >■ 7r"'"7r~ for CLEO 
II(II.V). We accept events with M within 5.2 - 5.3 GeV/c^ and \AE\ < 200(300) MeV for 
decay modes without (with) a 7r° in the final state. This fiducial region includes the signal 
region and a sideband for background determination. 

We have studied backgrounds from b c decays and other h ^ u and h ^ s decays and 
find that all are negligible for the analyses presented here. The main background arises from 
e^e~ — >• qq (where q = u,d,s,c). Such events typically exhibit a two-jet structure and can 
produce high momentum back-to-back tracks in the fiducial region. To reduce contamination 
from these events, we calculate the angle 6*5 between the sphericity axis of the candidate 
tracks and showers and the sphericity axis of the rest of the event. The distribution of 
cos ^5 is strongly peaked at ±1 for qq events and is nearly fiat for BB events. We require 
I cos ^5 1 < 0.8 which eliminates 83% of the background. Using a detailed GEANT-based 
Monte-Carlo simulation we determine overall detection efficiencies {S) of 11 — 45%, as 
listed in Table |. Efficiencies include the branching fractions for —>■ Kg tt+tt^ and 
7]-0 ryry -where applicable. 

Additional discrimination between signal and qq background is provided by a Fisher 



discriminant technique as described in detail in Ref . |]T0| . The Fisher discriminant is a linear 
combination, = X^i^i o^iVi, where the coefficients ctj are chosen to maximize the separation 
between the signal and background Monte Carlo samples. The 11 inputs, tji, are the cosine 
of the angle between the candidate sphericity axis and beam axis, the ratio of Fox- Wolfram 
moments H2/HQ [Q, and nine variables that measure the scalar sum of the momenta of 
tracks and showers from the rest of the event in nine angular bins, each of 10°, centered 
about the candidate's sphericity axis. 

We perform unbinned maximum-likelihood (ML) fits using AE, M, JF, the angle between 
the B meson momentum and beam axis, and dE/dx (where applicable) as input information 
for each candidate event to determine the signal yields. Four different fits are performed, 
one for each topology {h~^h~, h^ir^, h'^Kg, and K^n^, referring to a charged kaon or 
pion). In each of these fits, the likelihood of the event is parameterized by the sum of prob- 
abilities for all relevant signal and background hypotheses, with relative weights determined 
by maximizing the likelihood function (£). The probability of a particular hypothesis is 
calculated as a product of the probability density functions (PDFs) for each of the input 
variables. Further details about the likelihood fit can be found in Ref. []TU[. The parameters 
for the PDFs are determined from independent data and high-statistics Monte Carlo sam- 
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TABLE I. Experimental results. Branching fractions (B) and 90% C.L. upper limits are given 
in units of 10~^. The errors on B are statistical and systematics respectively. Reconstruction 
efficiency £ includes branching fractions of Kg — s- ■k'^tt~ and vr*^ — > 77. We quote the central 

value branching fraction in vr^vr'^ for convenience only. The statistical significance of the excess 
above background in this final state is insufficient for a first observation of this decay mode. 



Mode 


r* f fry \ 




Signif.(std.dev.) 


;B(10 ^) 




45 


A 7+1.8 
^- ' -1.5 


4.2 


4.7tJ|±0.6 




41 


5.4;2-i±1.5 


3.2 


< 12 


K+TT- 


45 


i8.8i^:^ 


11.7 


18.812:^ ± 1-3 




38 


12.ll|° 


6.1 


-19 -1+3.0+2.1 
-'-^•-'--2.8-1.4 




14 




7.6 


18.2l^-[; ± 1.6 




11 


14.811? 


4.7 




K+K- 


45 




0. 


< 2.0 




14 




1.1 


< 5.1 



pies. We estimate a systematic error on the fitted yield by varying the PDFs used in the 
fit within their uncertainties. These uncertainties are dominated by the limited statistics in 
the independent data samples we used to determine the PDFs. The systematic errors on the 
measured branching fractions are obtained by adding this fit systematic in quadrature with 
the systematic error on the efficiency. 

IV. RESULTS 

We summarize all branching fractions and upper limits in Table |. We find statistically 
significant signals for the decays B —>■ K'^n^, B tt+tt", B"^ K'^tc^, B"^ — > K^vr^, and 
B ^ KOtt". 

Fig. ^ shows the results of the likelihood fit for B -k'^-k'^ and B K'^tt^. The curves 
represent the na contours, which correspond to the increase in — 21n£ by n"^. The dashed 
curve marks the 3a contour; systematic uncertainties are not included in any contour plots. 
The statistical significance of a given signal yield is determined by repeating the fit with the 
signal yield fixed to be zero and recording the change in — 21n£. Fig. |^ shows distributions 
in M and AE for events after cuts on the Fisher discriminant and whichever of M and AE is 
not being plotted, plus an exclusive classification into i^vr-like and vrvr-like candidates based 
on the most probable assignment with dE/dx information. The likelihood fit, suitably scaled 
to account for the efficiencies of the additional cuts, is overlaid in the AE distributions to 
illustrate the separation between Ktt and vrvr events. 

We also compute from the PDFs the event-by-event probability to be signal or continuum 
background, and also the probability to be i^Tr-like or TTTr-like. From these we form likelihood 
ratios, 7^,,, = (P^, + P^J/(P;. + P^. + P^^ + P]^. + P^^^) and = P^J{P^^ + P^J. 
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FIG. 2. Contours of the -21n£ for the ML fit to Nk±t,^ and N^+^-, the K^ir"^ and tt+tt^ 
yields respectively. 
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FIG. 3. Projections of K-k and tttt events onto M and /S.E with cuts. Upper left: M distribution 
of KiT-Yike events; upper right: M distribution of TTTr-hke events. Lower left: /S.E distribution of 
events prior to tttt vs Ktt vs KK selection according to dE/dx; Lower right: /S.E distribution of 
events that are more likely to be tttt than Ktt or KK based on dE/dx. Overlays in the lower plots 
are the results of the likelihood fit scaled by the efficiency of the cuts used to project into these 
plots. Solid line: total fit; dashed: Kit; dotted: tttt; dot-dash: continuum background. 
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Superscript s and c stand for signal and continuum background respectively. Fig. ^ illustrates 
the distribution of events in TZsig (vertical axis) and 7^^ (horizontal axis). Signal events 
cluster near the top of the figure, and separate into Kir-like events on the left and vrvr-like 
events on the right. 

We find no evidence for the decay B — > K^K^ and set an upper limit accordingly, as 
shown in Table |I[ 

The results of the Kst^^ fit are shown in Fig. ^. The signal yield of IS.Slgg events is 
4.7(7 significant and robust under variations of cuts and PDF parameter variations. 

Figures ^ and |^ illustrate contour plots of —2 ln£ for the ML fit to K^h^ and h^n^. The 
branching ratios and limits associated with these four fits are given in Table |. We also show 
projections onto M and AE for these modes in Figures || and 

To evaluate how systematic uncertainties in the PDFs affect the statistical significance 
for modes where we claim first observations, we repeated the fits for the h'^h~ and K^n^ 
modes with all PDFs changed simultaneously within their uncertainties to maximally reduce 
the signal yield in the modes of interest. Under these extreme conditions, the significance 
of the first-observation modes tt+tt^ and K^tt^ becomes 3.2 and 3.8 a respectively. We 
also evaluate the branching ratios with alternative analyses using tighter and looser cuts 
on the continuum suppressing variable | cos6'5|. These variations correspond to halving and 
doubling the background in the fitted sample. The branching ratios change under these 
variations by much less than the statistical error. 

The fitted yields in the remaining modes are not statistically significant. We calculate 
90% confidence level (CL.) upper limit yields by integrating the likelihood function 



UL 



/O ^n.ax(iV)rfiV 

where Cma.x{N) is the maximum C at fixed to conservatively account for possible corre- 
lations among the free parameters in the fit. We then increase upper limit yields by their 
systematic errors and reduce detection efficiencies by their systematic errors to calculate 
branching fraction upper limits given in Table I. 



V. INFORMATION ON THE WEAK PHASE 7 

Charmless hadronic B decays are a sensitive probe of 7, the phase of the CKM matrix 
element V*ij, due to the interference of tree (Fig. |l|(a)) and penguin (Fig. 0(b)) diagrams. A 
number of methods for extracting 7 from these decays have been proposed, relying on either 
the construction of amplitude triangles 0], or ratios of CP-averaged branching fractions 

While some of these methods have the virtue of being independent of model assumptions 
about strong interaction effects none of them provide useful constraints on 7 given the 
present level of precision of our data, as they use only a restricted set of measurements. 

Alternatively, one may trade model independence against exhaustive use of the available 
data and attempt a model dependent fit to a large number of measurements. Such a fit is 
described in detail in Ref. In this section we will briefly describe the basic ideas, as well 
as the main results of this work. 
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' 0.00 0.25 0.50 0.75 1.00 



FIG. 4. The horizontal axis shows P^^/{P^^ + Pf^^) while the vertical axis depicts 

{P^^ + Pf^^)/{P^^ + Pf^^ + + P^^ + Pkk)- Superscript s and c stand for signal and con- 
tinuum background, respectively. Signal events cluster near the top of the figure, and separate into 
i^vr-like events on the left and TTTT-like events on the right. 
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B — > Ks7T° Likelihood Function 
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FIG. 5. i? — > KsTT^: Upper left: Likelihood function versus Kst^^ yield; Upper right: projection 
onto M axis (efficiency = 0.63); Lower left: projection onto AE axis (efficiency = 0.63); Lower 
right: projection onto J-" axis (efficiency = 0.69). Overlaid curves are the results of the likelihood 
fit, scaled according to the efficiency of the prcJj^ction cuts. Solid: total fit; dotted: continuum 
background; dashed: signal. 




FIG. 6. Contours of the -2 In/: for the ML fit to Nj^o^± and Nj^oj^±. 
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FIG. 7. Contours of the -21n£ for the ML fit to Nk±t,o and N^±^o. 
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FIG. 8. B —>■ KsTT^- Projections onto M (top) and AE (bottom) in GeV. Cuts are made on 
the Fisher discriminant and whichever of M and AE is not being plotted. 
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Mass for KttO candidates 



15 



n 



n 



n 





5.200 



5.225 



5.250 
Mass 



5.275 



5.300 



-0.30 



AE for KttO condidates 




-0.15 0.00 0.15 

AE kaon hypothesis 



0.30 



FIG. 9. B ^ K^TT^. Projections onto M (left) and AE (right). Cuts are made on the Fisher 
discriminant and whichever of M and AE is not being plotted. 



As pointed out in a recent paper by He, Hon, and Yang, [|T5l a number of CLEO measure 



ments of charmless hadronic B decays indicate a preferance for cos 7 < 0.0. This conclusion 
is based on a comparison of the CP-averaged branching fraction measurements with theo- 
retical predictions based on the factorization model. In this section we report results of a 
global fit that quantifies the qualitative observation made in Ref. [ISj . 



Assuming factorization one may express [|T^] [|T^ the decay amplitudes in terms of CKM 
matrix elements, form factors, decay constants, quark and meson masses, short distance 
coefficients, etc. Many of these are known to more than adequate precision. In addition, one 
may relate many of the poorly known quantities assuming unitarity of the CKM matrix as 
well as relationships among form factors. As a result, only five poorly known parameters are 
needed to predict CP averaged branching fractions for many of the CLEO measurements on 
charmless hadronic B decays. These parameters are: 

7 = Arg{V:,) 

\Vub/Vcb\ 

R = (2) 

{mt~mu)(ma+mu) m,i,xma ^ ' 

F^^'^ = B ^ 7C transitions form factor 
Aq^ = B ^ P transitions form factor 



We then form a between the CLEO results |T^ for the final states 



K^n^, fsT^vrO, K'^^tt^, K^n^, Tr+Tr", 7r±7r°, p^vr^, tU7r±, p^vr^, K*^7i^, ujK^^, coK^, ^K°, ^K^ 
and the theoretical predictions for the respective CP-averaged branching fractions. As addi- 
tional constraint we add iKife/Kbl = 0.08 ± 0.02, thus arriving at a 15 — 5 degree of freedom 
fit. We should note that only 9 of the 14 decay modes are unambiguously observed (i.e. 
with statistical significances of more than 4 standard deviations). For the remaining 5 final 
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states we see excess yields above background expectations that have statistical significances 
ranging from ~ 1.5 to close to 4 standard deviations. 

Our choice of decay modes to include was dictated by the following rationale. We consider 
for inclusion all decay modes with final states containing two pseudo-scalars, or a pseudo- 
scalar and a vector meson, for which we have results presented at this conference. We 
then exclude final states with 77 and 77' due to the known ambiguities in predicting these 
decays [jl8|. We furthermore exclude all final states for which the factorization model makes 
predictions that are well below the sensitivity of our present data. 

Minimizing the resulting we find a global minimum of 10.3 for 10 degrees of freedom. 
The parameters at the minimum are: 

7 = 113^1 (161) degrees 

IWKfel = 0.082 ±0.017 (0.062) 

Rsu = 1.7 ±0.4 (0.57) (3) 
pB^^ = 0.27l[!:[J^ (0.43) 
Aj^^' = 0.52t°oil (0.78) 

The central value for Rgu corresponds to ms{mh = 4.34GeV) of about 66 MeV. The 
numbers in parentheses are the values for a second local minimum with only marginally 
larger x^- These are the only minima found while repeating the fit 200 times with random 
starting points across a large volume in the five dimensional parameter space. 

Figure [1^ shows the dependence of versus 7. We refit at each value of 7 to allow the 
fit to find the global minimum at each fixed value of 7. This is important as the correlation 
matrix for the five free parameters has significant off-diagonal elements. For example, the 
two most correlated parameters are |Kjf,/V^b| and A^^^ with a correlation coefficient of —0.8. 

We note that the best fit values for all parameters are closely consistent with theoretical 
expectations |16| [|1^ |19|. This is either a surprising coincidence or an indication that the 



model provides an adequate description of nature. If we were to believe the latter than we 
would have to conclude that we have made a meaningful measurement of 7, the phase of 



VI. SUMMARY 

In summary, we have measured branching fractions for all four exclusive B Kir decays 
and made first observations of the decays B tt+tt^ and B K^n^. The latter observation 
completes the full set of i? — ^ Ktt measurements. In addition, we have shown that a global 
fit to CP-averaged branching fractions measured by CLEO allows for a model dependent 
measurement of 7 = II3I23 degrees. This is the first determination of the complex phase of 
the CKM matrix by any method other than the unitarity triangle construction. 
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FIG. 10. versus 7 curve for global fit to CP-averaged charmless hadronic B decay branching 
fractions. The triangles depict a local minimum that has the feature of = 218MeV. The squares 
depict the curve with the additional constraint of mg < 150MeV. 
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